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Abstract 

We prove existence of rotating star solutions which are steady-state solutions of the 
compressible isentropic Euler-Poisson (EP) equations in 3 spatial dimensions, with pre- 
scribed angular momentum and total mass. This problem can be formulated as a varia- 
' tional problem of finding a minimizer of an energy functional in a broader class of functions 

• having less symmetry than those functions considered in the classical Auchmuty-Beals 

CO . paper. We prove the nonlinear dynamical stability of these solutions with perturbations 

I having the same total mass and symmetry as the rotating star solution. We also prove 



' local in time stability of W ■°°(R'^) solutions where the perturbations are entropy- weak 



solutions of the EP equations. Finally, we give a uniform (in time) a-priori estimate for 
entropy-weak solutions of the EP equations. 
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1 Introduction 



The motion of a compressible isentropic perfect fluid with self-gravitation is modeled by the 
Euler-Poisson equations in three space dimensions (cf [1]): 

Pt + V- (pv) = 0, 

< (/9v)t + V • (pv V) + Vp(/9) = -/9V$, (1.1) 

A<I) = 47r/3. 

Here p, v = {vi,V2,v^), p{p) and $ denote the density, velocity, pressure and gravitational 
potential, respectively. The gravitational potential is given by 

$(x) = - / ,^^dy = -p*^ (1.2) 
Jr3 \x - y\ \x\ 

where * denotes convolution. The momentum pv is denoted by m = (mi, m2, w-s). System 
(jl.ip is used to model the evolution of a Newtonian gaseous star ( [4J ) . In the study of time- 
independent solutions of system (jl.ip . there are two important cases, non-rotating stars and 
rotating stars. A non-rotating star solution is a time-independent spherical symmetric solu- 
tion of the form (pat, 0, $Ar)(a;) (the velocity is zero), with <I>Ar(x) = —PN*-^y A rotating star 
solution models a star rotating around the xs-axis (x = (xi, X2, X3)) with prescribed angular 
momentum (per unit mass), or angular velocity. The existence and properties of stationary 
non-rotating star solutions is classical (cf. [4J). In contrast, the study for rotating stars is 
more challenging and of significance in both astrophysics and mathematics. A rigorous math- 
ematical theory for rotating stars of compressible fluids was initiated by Auchmuty &: Beals 
([1]) in 1971. The existence and properties of rotating star solutions were obtained by Auch- 
muty & Beals ([I]), Auchmuty ([2]), Caffarelli &: Friedman ([3]), Friedman &: Turkington([T3], 
[I3j), Li([2l]), Chanillo & Li([5|), and Luo & Smoller ([25]). In [26], McCann proved an 
existence result for rotating binary stars. 

The existence of rotating star solutions of compressible fluids was first obtained by Auch- 
muty & Beals (lU) who formulated this problem as a variational problem of finding a min- 
imizer of the energy functional F(p), (which will be defined in Section 2), in the class of 
functions Wm,s = Wm H Ws, where Wm is the set of integrable functions p -.M!^ —>■ which 
are a.e. non-negative, axi-symmetric, of total mass M = p{x)dx, and having a finite 
rotational kinetic energy (precise statements can be found in Section 2). Ws is defined by 

I^S = {p : ^ p(xi,X2,-X3) =p(xi,X2,X3), x^ G M, i = 1,2,3}. (1.3) 

In this paper, we first give a proof of the existence of a minimizer of the energy functional 
-F(p) in the wider class of functions Wm- Our proof is quite different from that in [1]. As in 
[1], the main difficulty in the proof is the loss of compactness due to the unboundedness of M^. 
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The method in [T] is to minimize the functional F on Wr = {p S Wm,s-,p{x) = > R} 
and to obtain some uniform estimates on the support of the minimizer. Our method is to use 
the concentration-compactness method due to P. L. Lions ([M])) which was also used in [29] 
to prove the existence of non-rotating star solutions. The reason that we seek minimizers in 
Wm instead of Wm,s is that we want to discuss the full stability problem dynamically in a 
more general context with less restrictions on the symmetry of solutions. 

The dynamical stability of these steady-state solutions is an important question. The 
linearized stability and instability for non-rotating stars and rotating stars were discussed by 
Lin ([23] ), Lebovitz ([12]) and Lebovitz & Lifschitz ([20]). The nonlinear dynamical stability 
of non-rotating star solutions was studied by Rein ([30j) via an energy-Casimir technique. 
It should mentioned here that the energy-Casimir technique was used in |16j to study the 
stability problem in stellar dynamics. Roughly speaking, for p[p) = , the result in [30] 
says that if the initial data of the Euler-Poisson equations is close to the non-rotating 
star solution in some topology, then the solution of (jl.ip with the same total mass as the 
non-rotating star, stays close to the non-rotating solution in the same topology as long as the 
solution preserves both the energy E[t) which is defined by 

E{t)=f (2^ + \p\^A{x,t)dx-^f \V^\x,t)dx, {I A) 

and the total mass p{x,t)dx. An interesting feature of the energy is that it has both 
positive and negative parts, making the analysis difficult. For solutions of (jl.ip without 
shock waves, energy is conserved. For solutions with shock waves, the energy E(t) is non- 
increasing due to the entropy condition associated with shock waves (cf. [18] and [32]). In 
this paper we extend the above nonlinear stability results to rotating stars. 

As in the non-rotating star case ([30]), our nonlinear stability result is in the class of solu- 
tions having the same total mass as that of the rotating steady-state solution. For solutions 
with different total masses, we investigate the nonlinear dynamical stability of a solution 
u = (/), v,^) G ^/o^) (which includes both rotating and non-rotating stars), in the context 
of weak entropy solutions, for more general perturbations not necessarily having the same 
mass as u, under some assumptions on the L°°-norm and the support of the solutions. This is 
achieved by using the techniques of relative entropies together with a careful analysis of the 
gravitational energy; i.e., the negative part in the total energy E{t). It should be mentioned 
here that the method of relative entropies was used by Dafermos ([9J) and Chen/Frid [6]) 
to study the stability and behavior of solutions of hyperbolic conservation laws. The main 
difficulty in applying this method to the the Euler-Poisson equations (jl.ip is again due to the 
non-definiteness of the energy density. We also give a uniform a priori estimate for the weak 
solutions of Cauchy problem of (jl.ip satisfying the entropy conditions. 

This paper is organized as follows: in Section 2, we prove the existence of rotating star 
solutions which are the minimizers of an energy functional F in Wm with prescribed total 
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mass and angular momentum with finite rotational kinetic energy. We also derive some 
properties concerning the minimizing sequence. These properties are interesting, and are 
important for our stability analysis. In Section 3, we prove our nonlinear stability result for 
rotating stars. Section 4 is devoted to the stability result for the entropy weak solutions and 
in Section 5, we obtain uniform in time a priori estimates for entropy weak solutions. 

Throughout this paper, for simplicity of presentation, we assume that the pressure function 
p{p) satisfies the usual 7-law, 

p{p) = P\ P> 0, (1.5) 

for some 7 > 1. We now introduce some notation which will be used throughout this paper. 
We use / to denote Jj^s, and use || • | |g to denote 1 1 • ||j;^q(]g3). For any point x = (xi, X2, X3) G M^, 
let 

r{x) = \J x\ + xl, z{x) = X3, Bnix) = {y e R^, \y - x\ < R}. (1.6) 
For any function / G L"'^(M^), we define the operator B by 

Bf{x)= [ JiyLdy = f*-^. (1.7) 

J \x-y\ \x\ 

Also, we use V to denote the spatial gradient, i.e., V = V^; = {dx^, dx^, dx^)- C will denote 
a generic positive constant. 



2 Existence of Rotating Star Solutions 



A rotating star solution (/j, v, <&)(r, z), where r = ^Jx\+ x'2 and z = X3, x = (xi, a;2, X3) G M^, 
is an axi- symmetric time-independent solution of system (jl.ip . which models a star rotating 
about the X3-axis. Suppose the angular momentum (per unit mass), J{mp{r)) is prescribed, 
where 

mp{r) = / p{x)dx = j 2tts j p{s, z)dsdz, (2-1) 

J ^/xf+x^<r Jo J -00 

is the mass in the cylinder {x = {xi,X2,x^) : \/x\ + x'2 < r}, and J is a given function. In 
this case, the velocity field v(x) = {vi,V2,V3) takes the form 

v(x) = (-"^^(7^"^\"^'^(7^"^\0). 
Substituting this in (jl.ip . we find that p{r,z) satisfies the following two equations: 

drP{p) = pdr{Bp) + pL{mp{r)r-^ , 

[2.2) 

dzP{p) = pdz{Bp), 
where the operator B is defined in (jl.7p . and 

L{mp) = J'^irup) 
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is the square of the angular momentum. For any function p > and 7 > 1, we define 

A(p) JL-. (2.3) 

7 — 1 7 — 1 

It is easy to verify that (cf. [T]) (j2.2p is equivalent to 

A!{p{x)) + / L{mp{s)s'^ds - Bp{x) = A, where p{x) > 0, (2.4) 

Jr{x) 

for some constant A. Here r(x) and z(x) are as in ()1.6p . In [Ij, Auchmuty and Beals formulated 
the problem of finding solutions of (j2.4p as the following variational problem. First, let M be 
a positive constant and let Wm be the set of functions p defined by (cf. (1.5)), 

Wm ={p : ^ M, p is axisymmetric, p>0, a.e., p G L'^{R^) n L'^(M^), 

/ pix)dx = M, / ' \ ^ dx < +00. j 
J J r{xy 

For p G Wm, we define the energy functional F by 

(^||5/)||2 < +00 follows from p € ^^(M^) n LT(M3) and Lemma 2.3 if 7 > 4/3.) In ([23]), the 
first term denotes the potential energy, the middle term denotes the rotational kinetic energy 
and the third term is the gravitational energy. Assume that the function L G C^[0,M] and 
satisfies 

L(0) = 0, L(m) > 0, for < m < M. (2.6) 

Auchmuty and Beals (cf. jT]) proved the existence of a minimizer of the functional F{p) in 
the class of functions Wm,s = Wm H Ws, where 

Ws = {p:M.^^ M, p{xi,X2, -X3) = p{xi,X2,X3), x,eR,i = l, 2, 3}. (2.7) 

Their result is given in the following theorem. 

Theorem 2.1. (111). If ^ > 4/3 and \2. 0\) holds, then there exists a function p{x) £ Wm,s 
which minimizes F[p) in Wm,s- Moreover, if 

G = {x(^R? : p{x) > 0}, (2.8) 

Then G is a compact set in M'^, and p G C^{G) Pi C^(M^) for some Q < j3 < 1. Furthermore, 
there exists a constant p < such that 

^'(Pi^)) + I^^) L{mp{s)s-^ds - Bp{x) = p, xeG, 
l^^^ L{mp{s)s~^ds - Bp{x) >p, X G M3 - G. 
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In this paper, we are interested in the minimizer of functional F in the larger class Wm- 
By the same argument as in [IJ, it is easy to prove the following theorem on the regularity of 
the minimizer. 

Theorem 2.2. Let p he a minimizer of the energy functional F in Wm and let 

T = {x£R^ : p{x) > 0}. (2.10) 

If ^ > 6/5, then p G C(M^) n C^{r). Moreover, there exists a constant A such that 

A'{p{x)) + l°^^^L{mp{s)s-^ds-Bp{x)=X, xGT, 
L{mp{s)s-^ds - Bp{x) > A, x € _ r. 

We call such a minimizer p a rotating star solution with total mass M and angular mo- 
mentum yLAjn). 

In this paper, we prove the existence of a minimizer for the functional F in the class Wm- 
For this purpose, in addition to (|2.6p . we require that L satisfies the following conditions: 

L{am) > a^/^L{m), 0<a<l, 0<m<M, (2.12) 

L'{m) > 0, < m < M. (2.13) 

Remark 1. Condition (j2.13p is called the Solberg stability criterion, see [33, Section 7.3]. This 
condition was also used by Auchmuty in ^ for the study of global branching of rotating star 
solutions. 

Our main result in this section is the following theorem. 

Theorem 2.3. Suppose that 7 > 4/3 and the square of the angular momentum L satisfies 
(KM), WdW and (KW- Then the following hold: 

(1) the functional F is hounded helow on Wm and infiy^^ P{p) < 0, 

(2) if {p^} C Wm is a minimizing sequence for the functional F, then there exist a sequence 
of vertical shifts a^es (ai G E, 63 = (0,0, l)j, a suhsequence of{p^}, (still labeled {p^}), and 
a function p G Wm, such that for any e > there exists R > with 

I p\x)dx > M - e, i E N, (2.14) 

Jaie3+Bfl(0) 

and 

Tp^(x) := p^(x + a^es) ^ p, weakly in L'^(M^), as i ^ 00. (2-15) 

Moreover (3) 

VB{Tp') VB{p) strongly in L'^{M.^), as i ^ 00. (2.16) 
(4) p is a minimizer of F in Wm ■ 
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Thus yO is a rotating star solution with total mass M and angular momentum yL. 

Remark 2. It is easy to verify that the functional F is invariant under any vertical shift, i.e., 
if p{-) G Wmi then p{x) =: p(x + aes) G Wm and F{p) = F{p) for any a G M. Therefore, if 
{/)*} is a minimizing sequence of F in Wa-/, then {Tp*} defined in (2.15) is also a minimizing 
sequence in Wm- 

Remark 3. In [13], [14] and [5], the diameter estimate of rotating star solutions with the 
symmetry p{r,—z) = p{r,z) was obtained. The ideas and techniques developed in [TH], [13] 
and [5] should also be applied to obtain the diameter estimates for the rotating star solutions 
in Theorem 12. 3i Due to the length of this paper, we leave this issue for the future study. 

Theorem 12.31 is proved in a sequence of lemmas. We first give some inequalities which will 
be used later. We begin with Young's inequality (see [l7j, p. 146.) 

Lemma 2.1. If f e r\U' , l<p<q<r< +oo, then 

ll/ll.<ll/llpll/ll^^ « = pl'"! - (2.17) 

The following two lemmas are proved in [Ij. 

Lemma 2.2. Suppose the function f G L^iR-^) n L'?(M3). If 1 < q < 3/2, then Bf =: f 
is in L^(M3) /or 3 < r < 3g/(3 - 2q), and 

115/11. <c(||/||5||/||^-^ + ||/||5||/||J-j, (2.18) 

for some constants C > 0, < b < 1, and 0<c<l. If q > 3/2, then Bf{x) is a hounded 
continuous function, and satisfies (2.18) with r = oo. 

Lemma 2.3. For any function f G L'^{R^) D L'^{R^), if j > 4/3, then VBf G L'^{R^). 
Moreover, 

I j f{x)Bf{x)dx\ = ^\\VBf\\l<c(^j \f\'"\x)dx^ (I \f\{x)dxY\ (2.19) 
for some constant C . 

Throughout this paper, we assume the function L, the square of the angular momentum 
satisfies conditions ^X^, ^XT^aud (2.13). Let 

/a/= inf F{p). (2.20) 

P&Wm 

We begin our analysis with the following lemma. 
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Lemma 2.4. Suppose 7 > 4/3. If p ^ Wm, then there exist two positive constants Ci and 
C2 depending only on 7 and M such that 



p{x)L{mp{r{x))) 
J V-;^ ^ 

This implies 



■p'{x) + ^^y±I±^!lEptlIl]dx < C,F{p) + C (2.21) 



fu > -00, 



where Jm is defined in (2.20). 

Proof. Using (j2.19p . we have, for p E Wm-, 



+ l p(':)L(^nMm ^^^ _ ci,2/3 ^1/3^^ (2.22) 
2 r(x)^ y 



Taking p=l,g = 4/3, r = 7, and a = 4^^^- in Young's inequality p.lTh . we obtain, 

\\p\\a/,<\\p\\1\\p\\\-^ = M'^\\p\\l^-\ (2.23) 

This is 

j p^/^dx < Mt"( J p^dx)\ (2.24) 
where h = s^^zriy- Since 7 > 4/3, we have < 6 < 1. Therefore, (|2.22p and (|2.24p imply 

\^iP) + < ^ _ i).Ml«+i( I A(p)d.)^ (2.25) 

Using ([2:21]) and the inequality (cf. [17J p. 145) 

al3 <ea' + e-^''pi\ (2.26) 

if s~^ + = 1 (s,t > 1) and e > 0, since 6 < 1, we can bound the last term in (j2.25p by 
^ j A{p)dx + C2, where C2 is a constant depending only on M and 7 (we can take e = 1/2 
and s = 1/6 and t = (1 — s~^)^^ in (j2.26p since s > 1 due to < 6 < 1). This implies 
(12:2111 . □ 

We also need the following lemma. 

Lemma 2.5. Suppose 7 > 4/3. T/ien 
(a) fu < /or every M > 0, 

(h) if [2JE) holds, then f^ > {M/Mf/^fu for every M > M >0 . 
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Proof. It follows from [T] that there exists p G Wm,s C Wm such that F{p) = ini p^Wm s ^(p)- 
By Theorem 12.11 it is easy to verify that the triple {p, v, $) is a time-independent solution 
of the Euler-Poisson equations (jl.ip in the region G = {x G : p{x) > 0}, where v = 

^_ x,J(rn,(r)) ^ a:,J(rn,ir)) ^ ^ Therefore 

^xP{p) = pV^{Bp) + pL{mp)r{x)-'^er: x G G, (2.27) 

where = (;^, 0). Moreover, it is proved in [3] that the boundary dG of G is smooth 
enough to apply the Gauss-Green formula (cf. [12]) on G. Applying the Gauss-Green formula 
on G and noting that p\qg = 0, we obtain, 

/ X ■ Va:Pip)dx = -3 / pip)dx = -3 / p{p)dx. (2.28) 
Jg Jg J 

By an argument in [33] (used also in [TO]), we obtain 

/ x ■ pV^Bpdx = -\ 1 pBpdx = -\ I pBpdx. (2.29) 
Jg ^ Jg ^ J 

(In fact, this can be verified as follows. Let 

1= [ X- pV.Bpdx = - [ p{x) I PMi^.Zyll^dydx. 
Jg Jg Jg \x-y\-^ 



Then 



/ = - / p(x) / PM-p^y^^dydx - f p{x) I P^^-^4^dydx 
Jg Jg \x-y\'^ Jg Jg F-yr 

p{x) [ _ , 

G Jg \x-y\'^ 



G 



pBpdx - /, (2.30) 
which is (j2.29p .) Next, since x • = r(x), we have 

X • p{x)L{mp{r{x))r^^ {x)erdx 

G 

p{x)L{m^{r{x))r^^ {x)dx 

G 

p{x)L{m^{r{x))r^'^{x)dx. (2-31) 



Therefore, from I^^B'UM^ we have 

— 3y p{p)dx = —- J pBpdx+ I p{x)L{mp{r{x))r^^{x)dx, (2.32) 

so that 

4 — 37 /■ 1 
F(p) = / p{p)dx / p(x)L(mp(r(x))r '^{x)dx. 

7 - 1 y 2 ' 
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Thus, if 7 > 4/3, F{p) < since L{m) > for < m < M. Since p £ Wm,s C Wm, then 
infp^y/M^ip) < 0- This completes the proof of part (a). 

The proof of part (b) follows from a scaling argument as in [29]. Taking h = (M/M)^/^ and 
letting p{x) = p{bx) for any p £ Wm- It is easy to verify that p € and that the following 
identities hold, 

pBpdx = / pBpdx, (2.33) 



Moreover, for r > 0, 



j A{p)dx = b-^ j A{p)dx. (2.34) 



oo 



mp{r) = 2-K s p{s,z)dsdz 
Jo J -oo 
pr poo 

= 27r / s p{bs, bz)dsdz 
Jo J ~oo 



'0 J -oo 

(•br foo 



2 rbr poo 

-r27r / s' p(s',z')ds'dz' 
b-^ Jo J-oo 



^mp{br). (2.35) 



Since L satisfies (j2.12p and 6 > 1, we have 



Thus, 



L{m-p{r)) > ^L{mp{br)). (2.36) 



p{x)L{m-p{r{x))) ^ 1 /■+°° 27rr , , ^^ , /, , n , , 
^ ^ ^ > ^ / —L{mp{br)) j p{br,bz)dzdr 

— oo 



r(x)^ 6^ Jq r 



-L{mp{r')) / p{r' , z')dz'dr' 



oo 



b^ Jq 

I m - 1 r( 'Till 

dx. (2.37) 



1 /■ p{x)L{mp{r{x))) 



b^ J r{xY 
Therefore, since 6 > 1, it follows from (|2.33p - (j2.37p that 

= {M/Mf/^F{p). (2.38) 
Since p ^ p is one-to-one between W^m and Wj^, this proves part (b). 



□ 



The following lemma gives the boundedness of a minimizing sequence of i*" in L'^( 
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Lemma 2.6. Suppose 7 > 4/3. Let {p*} C Wm be a minimizing sequence of F. Then {p*} 
is bounded in L'^(M^), and moreover, the rotating kinetic energy 

1 /■ p^(x)L(mp.(r(x))) 

2 J r(x)2 

is a/so uniformly bounded. 
Proof. By Lemma 12.41 we have 

The lemma follows from this and Part a) in Lemma 12.51 □ 

Lemma 2.7. Suppose 7 > 4/3. iyei {p*} C Wm be a minimizing sequence for F . Then there 
exist constants rg > 0, (5o > 0, io ^ N and G M'^ to^/i r(2;*) < r^, such that 

p\x)dx > 5o, i> io- (2.40) 

Bi(x^) 

Proof. First, since limj^oo -^(p*) — > /m and /a/ < (see part (a) of Lemma [23]) . for large i, 

- ^ < -F{p') p'Bp^dx. (2.41) 

For any i, let 

Now 

p''Bp\x)dx 



sup / p\y)dy. (2.42) 

xeK3 ./|y-x|<i 



p [x) I rdyax 

J«i \y - x\ 

p\x) [ -^^^dydx+ [ p\x) [ -^^^dydx+ [ p\x) [ -^^^dydx 

' J\y-x\<l \y - X\ Jr» Jl<\y-x\<r \y - x\ Jr^ J\y-x\>r IV - x\ 

=:Bi + B2 + B3, (2.43) 

and B3 < M^r-^. The shell I < \y — x\ < r can be covered by at most Cr^ balls of radius 1, 
so B2 < CMSir^ . By using Holder's inequality and (|2.18p . we get 

Bi = / p\x) / I I dydx 

Jr3 7r3 \y - x\ 

< \\p%/3\\B{x{y\\y-x\<i}{y)p\y)}\U 

< c\\p1,/s (\\x{yiiy-x\<i}{y)p\y)\\>%Js' + \\x{y\^^^^^^ 



< 



C\\p%/3 {S^\\P%J,' + mPllj^) , (2.44) 



11 



where x is the indicator function, < & < 1 and < c < 1. By lemma 2.6, we know that 
11/3*11^ is bounded, so ||/0*||4/3 is bounded if 7 > 4/3 in view of (2.17) and the fact ||p*||i = M. 
This gives Bi < C {5^ + 51). It follows that we could choose r so large that the above estimates 
give / p'^Bp^{x)dx < —Jm if were small enough. This would contradict (j2.4ip . So there 
exists (^0 > such that 6i > 60 for large i. Thus, as i is large, there exists G and io S N 
such that 

/ p\x)dx >do, i> io. (2.45) 

We now prove that there exists ro > independent of i such that those must satisfy 
r(x*) < ro for i large. Namely, since has mass at least 5o in the unit ball centered at x*, 
and is axially symmetric, it has mass > Cr(x*)5o in the torus obtained by revolving this ball 
around xs-axis (or z-axis). Therefore r{x'^) < {C6o)~^M. □ 

In order to prove Theorem 12.31 we will need the following lemma which is proved in [29] , 
and uses a concentration-compactness argument. 

Lemma 2.8. Suppose 7 > 4/3. Let {/*} be a bounded sequence in L'^(W^) and suppose 

f f weakly in L^{W^). 

Then 

(a) For any R> 0, 

yB{xBn{0)n ^ Vi3(xB«(0)/°) Strongly in L\W^), 
where x is the indicator function. 

(b) If in addition {/*} is bounded in L'^i^), f^ G L-'^(M'^), and for any e > there exist 
R > and G N such that 

[ \f\x)\dx<e, i>io, (2.46) 

J\x\>R 

then 

VBf' VBf^ strongly in L^i'M?). 

Before giving the proof of Theorem 12.31 we first outline the main steps. In step 1, we first 
show (2.15) and (2.16). In step 2 we show that if p is a weak limit in L'''(]R^) of {Tp*}, then 
mp{r) is a continuous function of r for all r > 0. The third step is to prove that F is lower 
semi-continuous with respect to the weak topology in L'^(M^). 

Proof of Theorem \2.3\ 

Step 1. We prove (2.16), and apply Lemma [THl to prove (2.14). We begin with a splitting as 
in [29]. For p £ Wm, for any < i?i < R2, we have 

P = PX\x\<Ri + PXRi<\x\<R2 + PX\x\>R2 =■ Pi + P2 + P3, (2.47) 
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where x is the indicator function. It is easy to verify that 

3 



j A{p)dx = Y,j A{p,)dx, (2.48) 



and 



where 



Also, 



J pBpdx = Y.[ PjBpjdx + /i2 + /i3 + I23, (2.49) 



p{x)L{mp{r{x)) f Pj{x)L{mp^{r{x))_^^ 



r2(x) ~^ J r2(x) 

(r 

r2(a;) 



f Pi(a:)(^(mp(r(x)) - L(mp^ (r(x)) 



It follows from (IXiTD - ilCTD that 

F(p)=X;F(p,)-i J] 

j=l l<«<j<3 

Since p > pj , we have mp{r) > nip. (r) for any r > and j = 1, 2, 3. By ()2.13p . 

F(p)>^i^(Pi)-- (2.52) 

j=l l<i<j<3 

Using (j2.52p and Lemma [231 by the same argument as in the proof of Theorem 3.1 in [29], 
we can show that 

/m - F{p) < C/mMiMs + C{R^' + ||p||^'?+i)/6||Vi?p2||2), (2.53) 

by choosing R2 > 2Ri in the splitting ()2.47p . where Mi = / pi{x)dx = Jj^|^^^ p(x)dx, 
M3 = / p-i{x)dx = p{x)dx and (7 = 1/(7 — 1). Let {p*} be a minimizing sequence of F 

in Wm- By Lemma |2.7| we know that there exists G N and 6q > independent of i such 
that 

/ p\x)dx >6o, i>io (2.54) 
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for some E with r(x*) < tq for some constant tq > independent of i. Let Oj = z(x*) 
and i?o = ro + 1, then (j2.54p imphes 

/ p'{x)dx > 60, if i > io, (2.55) 

where = (0,0,1). Having proved (j2.55p . we can follow the argument in the proof of 
Theorem 3.1 in ^ to verify U^Mti for 

rix)=Tp\x) =:p\. + a,es) 

by using (j2.52p and (j2.55p and choosing suitable Ri and R2 in the splitting (j2.47p . We sketch 
this as follows. The sequence Tp* =: /?*(• + a^es), i > io, is a minimizing sequence of F in 
Wm (see Remark 2 after Theorem 12. 3p . We rewrite (|2.55p as 

/ Tp\x)dx > (5o, i> (2.56) 
J Bum 

Applying (j2.53p with T/o* replacing p, and noticing that {Tp*} is bounded in L'''(]R'^) (see 
Lemma 2.6), we obtain, if R2 > 2Ri, 

- CfMMiMl < C{R^^ + llVSrp^lla) + F{Tp') - /m, (2.57) 

where M{ = jTp\{x)dx = j^^^^^J p\x)dx„ M| = jTp\{x)dx = j^^^^^Jp%x)dx and 
Tp2 = XRi<\x\<R2'^P^- Since {Tp"^} is bounded in L'''(M^), there exists a subsequence, still 
labeled by {Tp*}, and a function p £ Wm such that 

Tp' p weakly in L^{W^). 

This proves (2.15). By (p36]) . we know that M\ in (f237j) satisfies M\ > 5q for i > io by 
choosing Ri > Rq where Rq is the constant in ()2.56p . Therefore, by (j2.57p and the fact that 
/m < (cf. Part (a) in Lemma 2.5) , we have 

- C/mSoMI < CR^^ + C\\VBp2\\2 + C\\VBTpi - V^palb) + F{Tp') - /m, (2.58) 

where p2 = X\x\>R2P- Given any e > 0, by the same argument as [29], we can increase 
Ri > Rq such that the second term on the right hand side of (j2.58p is small, say less than 
e/4. Next choose R2 > 2Ri such that the first term is small. Now that Ri and R2 are fixed, 
the third term on the right hand side of (|2.58p converges to zero by Lemma 12.8( a). Since 
{Tp*} is a minimizing sequence of F in Wm, we can make F{Tp'^) — /m small by taking i 
large. Therefore, for i sufficiently large, we can make 

Ml =: [ Tp\x)dx < e. (2.59) 
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This verifies (j2.46p in Lemma 2.8 for /* = Tp*. By weak convergence we have that for any 
e > there exists R > such that 

M - e < / p{x)dx < M, 
Jbr{0) 

which imphes p £ L^(R^) with J pdx = M. Therefore, by Lemma l2.8r b) .we have 

\\VBTp' - VBp\ I2 ^ 0, +00. (2.60) 

This proves (2.16). (2.14) in Theorem O fohows from ([23^]) by taking R = R2. 

Step 2. Let p be a weak hmit of a subsequence of {Tp^} in L'^(IR^) (we stih label the subse- 
quence by {Tp*}). We claim that the mass function 

mp{r) =: / p{x)dx is continuous for r > 0. (2-61) 

This is proved as follows. By the lower semicontinuity of norms (cf. |22] p.51) and Lemma 
2.6, we have 

IIpII^ < lim \id\\Tp% = lim inf < C, (2.62) 

i— >oo i^oo 

for some positive constant C. For any e > 0, by the weak convergence and ()2.14p which we 
have already proved, there exists R > such that 

/ Tp\x)dx < e, i G N, (2.63) 

J\x\>R 

I p{x)dx = lim / Tp\x)dx < e. (2.64) 

J\x\>R «^°°ibl>fl 

For any r > and ri > r, 

< rnp{ri) — mp{r) 

p{x)dx 



and 

'\x\>R ^^'^J\x\>R 



/ p(x)dx + / p{x)dx. (2.65) 

J r<^J x\+x\<rx,\xz\>R J r<^fx\+x\<Tx,\xz\<R 

Since \x = (xi, X2, X3) G : r < \/x\+li^ < ri, {x^l > R} C {x = (xi, 2:2, ^3) G : \x\ > 
R}, by (IZMD . we have 

p{x)dx < e. (2.66) 



r<^/x'f+xl<ri,\xs\>R 

By (I2.62P and Holder's inequality, 



/ p{x)dx 

J r<^ xl+xl<ri,\xz\<R 

I \ 

{x = (xi,X2,X3) £W^:r<^Jxl+xl<ri, \x^\ < R} \ 



< WPW-y I meas 



< C[2TrR{ri + r)(ri - r)]^/^', (2.67) 
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where meas denotes the Lebsgue measure and 7' = (7— 1)/7. Now, if we take 6 = min{ 2-KR{2r+i) - 
then by ()2.67p . we obtain 

p{x)dx < e, (2.68) 



r<y/x'f+xl<ri,\x3\<R 

whenever < ri — r < (5. It follows from (j2.65p . (j2.66p and (j2.67p . we have 

\mp{ri) - mp{r)\ < 2e, (2.69) 

whenever < ri — r < 6. This proves that mp{r) is continuous from the right for any r > 0. 
By the same method, we can show that mp{r) is continuous from the left for any r > . 
Since mp(0) = 0, this proves ()2.6ip . 

Step 3. Let {p*} be a minimizing sequence of the energy functional F, and let p be a weak 
limit of {Tp*} in L'''(]R^). We will prove that p is a minimizer of F in Wm'-, that is 

F{p) < lim inf F{Tp'). (2.70) 

i— >oo 

First, by (|2.62p . we have 

/ A{p)dx < lim inf / A{Tp')dx. (2.71) 

J J-^OO J 

We fix a positive number 5 and show that 

r Tp^(x)L(mm^i(r(x)) — p(x)L(rn^(r(x)) 

/ ^ ^ ^ ^ ^ \;; , ^ dx = 0. 2.72 

•Jr(x)>s r-^ix) 

To see this, we write 

r {Tp'{x)L{'mJ'p^ (r(x)) - p{x)L{mp{r{x)) 
Jr{x)>5 r'^ix) 

I {Tp'jx) - p{x))L{mp{r{x)) 
Jr{x)>& r'^{x) 

Tp^{x){L{mrppi{r{x)) — L{mp{r{x))) 
r^{x) 

{Tp^{x) — p{x))L{mp{r{x)) 



lim 

i^oo Jr(x)>6 



+ / ^ /; ^ ' '" dx. (2.73) 

Jrix)>S r^[x) 



For any i? > 0, we have 

f (T^n^fT.) — n( Tn\^ T,( m.-( r( TfYi 

-dx 



r{x)>& r'^ix) 

(Tp'ix) - p{x))L{mp{r{x)) 

r{x)>5,\x\<R \X) 

{Tp'jx) - p{x))L{mp{r{x)) 

r(x)>S,\x\>R r2(x) 

In view of ()2.63p and (|2.64p . for any e > 0, we can choose R such that 
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/ 

Jr(: 



{Tp\x) - ~p{x))L{mp{r{x)) ^ 2L{M)e ^ ^^.75) 

lr{x)>&,\x\>R r'^ix) ~ 6^ ' 

By the weak convergence of {Tp^} in L'>'(R^) and the fact that L is defined on a bounded range, 
L{mp{r{x))x{r{x)>5,\x\<R}ix)'''^'^ix) ^ L'^'(M^), where as before x is the indicator function, 
and 7' = (satisfying 1/7 + 1/7' = !). We have 

^.^ f {Tp\x) - p{x))L{mp{r{x)) 

Jr{x)>5,\x\<R r^ix) 

= 0, (2.76) 
because Tp* converges weakly to p. Since e is arbitrary, (|2.75p and (12.761) imply 



r{x)>5 

We handle the second term in (j2.73p as follows. By weak convergence, we know that mrppi{r) 
converges to mp{r) pointwise for r > 0. Since mj>pi(r) and mp{r) are non-decreasing functions 
of r for r > and mp[r) is continuous on [0, +00) (see (I2.6ip ). by a variation on Dini's theorem 
(|31j. p. 167)*, we know that mrppi{r) converges to mp{r) uniformly on the interval [0,R] for 
any R> 0. Since L £ C^[0,M], it follows that L(mj'pi(r)) converges to L{inp{r)) uniformly 
on any interval [0, R]. For any e > 0, we can fix i? > such that (|2.63|) and (|2.64|) hold. Since 
L(mrppi{r)) converges uniformy to L{mp{r)) on any interval [0, i?], we have 

hm \\L{mTp^{-)) - L{m~p{-))\\L^[m = (2-78) 

Let 

^5 = {x G R^r(x) > 5}, (2.79) 

then we have, using (|2.63p and (I2.64p that 

Tp''{x){L{mTp^{r{x)) - L{mp{r{x))) 
r{x)>5 r\x) 

f Tp'{x){L{mT,{r{x))-L{mp{r{x))) 

< / 2^^^ dx\ 

f Tp'{x){L{mTpi{r{x)) - L{mp{r{x))) 
+ 1 / ^ dx\ 

< \\L{mTp.{-)) - L(mp(-))||L-[o,i?]'^~'M + 26-'L{M)e. (2.80) 



* We thank Dmitry Khanvinson for pointing out this to us. 
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Since e is arbitrary, it follows from (j2.78p and (j2.80p that 

T p''{x){L{m'rpi{r{x)) - L{mp{r{x))) 



2— »00 



lim / o/!. dx = 0. (2.81) 



This, together with ([2:73]) and (pTTl) . implies (p:72|) . Next, we show that 

' ^ ^ \;; , '^^ ^ ^' dx > 0, (2.82) 

r^[xj 



t—*oo 



by using (j2.72p and the monotone convergence theorem for integrals. In fact we have 

-dx 



T p'^{x)L{mxpx{r{x)) — p{x)L{mp{r{x))) 



r^{x) 

dx 



Tp'{x)L{mTpr{r{x))){l - XA, 



r'^{x) 

[Tp\x)L{mTpi{r{x))) - p{x)L{mp{r{x)))]xAi, , 

+ I n^^ dx 

r^[x) 

+ / P(^)^(^P(K^)))(XA.-1) ^^^ 



where x is the indicator function, and As is the set defined in (|2.79p . For any i > 1, 

-dx > 0. (2.84) 



Tp'{x)L{mTpi{r{x))){l - X Asj 
r'^[x) 

We fix 6, and by (|2.72p . we know that the second term on the right hand side of (|2.83p 
approaches zero as i — > oo. Therefore, in view of (|2.84p . 

, Tp'{x)L{mT .{r{x))) - p{x)L{mp{r{x))) 
hm mt / TTT—. dx 



r'^{x) 



> I Kx)L{mMxmXAs-l) ^^^ (2.85) 



By the monotone convergence theorem of integrals, we have 

lim I / Kx)HrnMxmXA,-l) ^ ^ 
5^0 J r^{x) 

Letting ^ in I^Ml, gives (12:821) . By ([M]) . (2.71) and (12:821) . we obtain 

F{p) < lim inf F{Tp'). (2.87) 

i—*oo 

Since Tp* is a minimizing sequence, p is a minimizer of F in Wm- This completes the proof 
of Theorem [ 
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3 Nonlinear Stability of Rotating Star Solutions 



We consider the Cauchy problem for (jl.ip with the initial data 

p{x, 0) = po{x), v(x, 0) = vo(x). (3.1) 
We begin by giving the definition of a weak solution. 

Definition: Let pv = m. The triple (p,m,$)(x,t) {x G R^,t £ [0,T]) (T > 0) and $ 
given by <\1.2\ with p> 0, m, m ® m/p and p\7^ being in Lj^^{M!^ x [0,T]), is called a weak 
solution of the Cauchy problem (jl.ip and (j3.ip on M'^ x [0, T] if for any Lipschitz continuous 
test functions and ^ = (V'l; "02 5 "^s) with compact supports in M3 X [o,r], 

J (pi/jt + m • V^) dxdt + J po{x)'il^{x, 0)dx = 0, (3.2) 

and 

^ y (^m • *t + "^'^"^ • V*^ dxdt + y" mo(x)*(x,0)(ix = J P^^^dxdt, (3.3) 
both hold. 

For any weak solution, it is easy to verify that the total mass is conserved by using a 
generalized divergence theorem for U functions (r > 1) (cf. [7]), 

/p(.,t)<i. = /.(..0)<i., *>0. (3.4) 

The total energy of system (jl.ip at time t is 

E{t) = E{p{t),v{t)) = I (^A{p) + Ip\v\^) ix,t)dx [ \V<i>Wx,t)dx, (3.5) 

where as before, 



2"' ' y Svr 



MP) = (3.6) 

7-1 

Note that the energy E(t) has both a positive and a negative part. This makes the stability 
analysis highly nontrivial, as noted in [30]. For a solution of (jl.ip without shock waves, the 
total energy is conserved, i.e., E(t) = E(0) {t > 0)(cf. [33]). For solutions with shock waves, 
the energy should be non-increasing in time, so that for all t > 0, 

E{t) < E{0), (3.7) 

due to the entropy conditions, which are motivated by the second law of thermodynamics (cf. 
|18j and |32j). This will be proved in Theorem 5.1, below. 
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We consider axi-symmetric initial data, which takes the form 
po{x) = p{r,z), 



^oix) = Vo(r, z)er + Uo(r, z)eg + Vq{p, z)e3. (3i 



Here r = A/xf + aJf, z = X3, x = {xi,X2,X3) G (as before), and 

e, = (xi/r, X2/r, Of, eg = (-xa/r, x^/r, Of, eg = (0, 0, 1)^. (3.9) 
We seek axi-symmetric solutions of the form 
p{x,t) = p{r,z,t), 

v(x, t) = v''{r, z, t)er + v^{r, z, t)eg + v^{r, z, t)e3, (3.10) 

$(x, t) = ^>(r, z, t) = -Bp{r, z, t), (3.11) 

We call a vector field u(x,t) = {ui,U2,U3){x) (x G M'^ ) axi-symmetric if it can be written in 
the form 

u(x) = u^{r, z)er + u^{r, z)eg + u^{p, 21)63. 

For the velocity field v = {vi,V2,V3){x, t), we define the angular momentum j{x, t) about the 
X3-axis at (x, t) , t > 0, by 

j{x,t) = X1V2 - X2V1. (3.12) 

For an axi-symmetric velocity field 

v(x, t) = v'^{r, z, t)er + v^{r, z, t)eg + v^{p, z, t)e3, (3.13) 

Vl = V V , V2 = V H V ,V3 = V , (3-14) 

so that 

j{x,t) = rv'^{r,z,t). (3.15) 
In view of ( I3.13P and ()3.15p , we have 

Therefore, the total energy at time t can be written as 

E{p{t)Mt)) = I A{p){x,t)dx + ^ J ^^^dx 

- — [ \VBp\'^{x,t)dx + - [ p{\v''\'^ + \v^\'^){x,t)dx. (3.17) 
8tt J 2 J 

There are two important conserved quantities for the Euler-Poisson equations (jl.ip : namely 
the total mass and the angular momentum. In order to describe these, we define Dt, the non- 
vacuum region at time t > of the solution by 

Dt = {x eR^ : p{x, t) > 0}. (3.18) 
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We will make the following physically reasonable assumptions Al)-A4) on weak solutions of 
the Cauchy problem (|l.ip and (j3.ip : 

Al) For any t > 0, there exists a measurable subset Gt C Dt with meas{Dt — Gt) = {meas 
denotes the Lebsegue measure) such that, for any x £ Gt, there exists a unique (backwards) 
particle path ^(r, x, t) for < r < t satisfying 

dri{T, X, t) = v(e(T, X, t), r), X, t) = X. (3.19) 

For X G Gt, we write 

e(o,x,i) = e-t(x). 

Also, for X G M'^ and t > 0, we denote the total mass at time t in the cylinder {y G : 
r{y) < r(x)} by mp(t)(r(x)), i.e., 



rripitMx)) = / p{y,t)dy. (3.20) 

For axi-symmetric motion, we assume 
A2) 

mpit){r{x))=mp,{r{^^t{x))), foTxeGt,t>0. (3.21) 
Moreover, the angular momentum is conserved along the particle path: 
A3) 

j(x,t) =i(C_t(x),0), iovx£Gt,t>0. (3.22) 
(Both (3.21) and (3.22) are shown in [33j if the solution has some regularity.) 

Finally, for L = f' , we need a technical assumption; namely, 
A4) 

L(mo(-f\(r) +mrt(r))mo-rti(r) 
lim ^ ^^^^ ' " ^ '^ ' = 0, (3.23) 

for t > 0, where a{t) = p{t) — p. 

Remark 4. (j3.23p can be understood as follows. For any p G Wm, we have limr^o+ rnpi'i') = 0. 
Therefore linv_»o+ -^("T'p(t)('^) + iT^pi^)) = L{0) = 0, so if we define 

r+oo 

p{s,t) -~p{s) = / {p{s,z,t) - p{s,z))dz, 



£ L°°(0, <5) /or some 6 > 0, (3.24) 



then if 

'm-a{t){r) _ jQ{2TTs{p{s,t) - ^{s))ds 



(|H:25]) will hold. If p{-,t) ->(•) G L°°(0, 5), then holds. This can be assured by assuming 

that p{r, z, t) — p{r, z) G L°^((0, 5) x M x M+) and decays fast enough in the z direction. For 
example, when p{x,t) — p{x) has compact support in and p{-,t) — p{-) G L°°(]R'^), then 
([3:23]) holds. 



21 



Now we make some assumptions on the initial data; namely, we assume that the initial 
data is such that the initial total mass and angular momentum are the same as those of the 
rotating star solution (those two quantities are conserved quantities). Therefore, we require 

Ii) 

J poix)dx = J p{x)dx = M. (3.25) 

Moreover we assume 

I2) For the initial angular momentum j(x, 0) = rvQ{r, z) =: jo{r, z) (r = \/ x'f + x'2, z = x^ 

for X = {xi,X2,X3), we assume j{x,0) only depends on the total mass in the cylinder {y G 
^^,r{y) < r{x)}, i.e. , 

jix,0)=joimp,irix))). (3.26) 

Finally, we assume that the initial profile of the angular momentum per unit mass is the same 
as that of the rotating star solution, i. e., 

I3) 

jo(m)=L(m), 0<m<M, (3.27) 

where L{m) is the profile of the square of the angular momentum of the rotating star defined 
in Section 2. ( (j3.26p implies that we require that VQ{r, z) only depends on r.) 
In order to state our stability result, we need some notation. Let A be the number in Theorem 
2.2, i.e., 

A'{p{x)) + L{mp{s))s-^ds - Bp{x) = A, x G F, 
l^^^ L{mp)is))s-^ds - Bp{x) > A, x G - F, 
with A defined in (2.3), and F defined in (2.10). 



(3.28) 



For peL^r\ LT, we define, 

d{p, p)= [ [A{p) - A{p)] + {p-p) r { ^^"^'j'^^ ds - A - Bp}dx. (3.29) 



' r{x) 

Remark 5. For x G F, in view of (j3.6p and (3.28), we have 

L(mp(s)) 



{A{p) - A{p)){x) + ( / -^^^ds - A - Bp{x)){p - p) 

Jr{x) ^ 

= (A{p)-A{p)-A'{p){p-p)){x) 

^ P{p)-P{~p)-P'mp-~P) ^^^^^_ (3.30) 
7 - 1 

Thus, for p G Wm^ 

d{p,p)>0. (3.31) 
Moreover, d{p, p) = if and only if p = p, and if 7 < 2, 

d{p,p)>C\\p-p\\l peWM. (3.32) 
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We also define 

, . 1 f p{x)L{mp{r{x)) - p{x)L{mii{r{x)) 
di{p,p) = 2 / dx 

L{mp{s))ds{p{x) — p{x))dx, (3.33) 

r{x) 

for p G Wm- We shall show later that di > 0. Our main stability result in this paper is the 
following global-in-time stability theorem. 

Theorem 3.1. Let p be a minimizer of the functional F in Wm, cin-d assume that it is unique 
up to a vertical shift. Suppose 7 > 4/3 and the above assumptions Al)-A4) and Ii)- I3) hold. 
Moreover, assume that the angular momentum of the rotating star solution p satisfies i2. 6\) . 
112.1^) and (2.13). Let {p,v,^){x,t) be an axi- symmetric weak solution of the Cauchy problem 
nn}) . lO) with p{-,t) e L^n L-y, p|v|2(.,t) E L^ and V$(-,t) = -VBp{-,t) e L^. If the 
total energy E{t) (c.f. (3.5)) is non-increasing with respect to t, then for every e > 0, there 
exists a number (5 > such that if 

diPo,p) + ^W'^Bpo - VBpWl + \di{po,p)\ 

+ \j Po{x){\vl\'' + \vl\''){x)dx < 5, (3.34) 
then there is a vertical shift ae^ (a G M, es = (0, 0, 1) ) such that, for every t > 
d{p(t),T'^p) + ^\WBp{t) - BT-PWI + |di(p(t),rV)| 
+ ^ y p{x, t){\v''{x, i)P + W{x, t)\^)dx < e, (3.35) 
where T"'p{x) =: p{x + aea). 

Remark 6. The vertical shift 063 appearing in the theorem is analogous to a similar phe- 
nomenon which appears in the study of stability of viscous traveling waves in conservation 
laws, whereby convergence is to a "shift" of the original traveling wave. 

Remark 7. Without the uniqueness assumption for the minimizer of F in Wm, we can have 
the following type of stability result, as observed in [30] for the non-rotating star solutions. 
Suppose the assumptions in Theorem 13.11 hold. Let Sm be the set of all minimizers of F in 
Wm and (p, v, $)(x, t) be an axi-symmetric weak solution of the Cauchy problem (jl.ip . ()3.ip 
with p{-,t) E L^nLT, />|v|2(.,t) G L^ and let V$(-,t) = -VBp{-,t) £ L^. If the total energy 
E{t) is non-increasing with respect to t, then for every e > 0, there exists a number 6 > 
such that if 

1 



inf 



dipoJ) + ^W^^Po - ^Bp\\i + \di{po,p)\ 



+ I I Po(x)(|^5l' + \yl\^)ix)dx < 6, (3.36) 
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then for every t > 

1 



inf 



d{p{t),T'^p) + ^llVSp(t) - BT^pWi + \d,ipit),T^p)\ 



+ \j t){\v\x, t)\^ + \v\x, t)\''){x)dx < e. (3.37) 

The proof of this follows exactly along the same line as that for Theorem 13.11 
In order to prove Theorem 13. 11 we need several lemmas. First we have 

Lemma 3.1. Suppose the angular momentum of the rotating star solutions satisfies i2. 6]) , 
(EdM) and (2.13). For any p{x) G Wm, if 

lim L{mp{r) + mp{r))mcj{r)r~'^ = (3.38) 

where a = p — p., then 

diip,p)>0, (3.39) 



where di is defined by < fX 



Proof. First, we introduce some notation. For an axi-symmetric function f{x) = f{r,z) 
(r = \/x\ + x'2, z = X3 for x = (xi, j;2, 2:3)), we let 

/+00 
f{r,z)dz, (3.40) 
-00 



mf{r) = I f{x)dx = / f{s)ds, 

J {x:yjx\+xl<r} Jo 



(3.41) 



so that 

m'f{r) = f{r). (3.42) 

In order to show (j3.39|) . we let 

a{x) = {p-p){x), (3.43) 

and for < a < 1, we define 

1 f {p + aa){x)L{mp+aair{x))) - p{x)L{mp{r{x))) 



Q(») = 2 J ^) 



00 

3; 



Then 



Since 



a I s ■^L{mp{s))dsa{x)dx. (3.44) 

' r(x) 



Q(0) = 0, Q{l)=di{p, p). (3.45) 

fr{x) f+00 

mpj^a(j{r{x)) = 1 2tts i {p + aa){s, z)dzds, (3.46) 

Jo J -00 
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we have 

rr{x) r+oo 



^ rr(x) f-+co 

—mpj^cicj{f{x))= / 2tts / a{s,z)dzds = mtj{r{x)). (3.47) 

"Ck Jo J-oo 



Therefore, 



Q'{a) = i / -i^)Lirnp-^^.{r{x))) 
2 J r (x) 

1 /■ {p + aa){x)L' {mp+aa{r{x)))ma{r{x)) 
+ 2 J r2(x) 



and in view of (j3.42p . 



s-'^L{mp{s))dsa{x)dx, (3.48) 



^L{mp+au{r)) = L'{mp+aair)){ p + aa){r). (3.49) 



Therefore, by virtue of (j3.49p and (j3.42p . we obtain 



1 f [p + a(j){x)L' {'mpj^aa{r{x)))m„{r{x)) 



+ 00 

( /5 + a(T)(r)L'(mp+a^(r))m^(r)r"^(ir 





2 7 r2(x) 
_ 1 

1 /■+°° d 



-^[L(mp+oCT(r))]m^(r)r ^dr. (3.50) 

For < a < 1, since (cf. (2.13)) L'{m) > 0, we have 

L{mp+aa{r)) < L{mp+p{r)). (3.51) 
This, together with (I3.38p . imphes 

lim L{mpj^acr{r))mcr{r)r~'^ = 0. (3.52) 

Moreover, since mo-(+cxD) = J a{x)dx = Jip — p){x) = and 

hm Lirup+aair) = L{M), 

we have 

lim L{mp+aa{r))m„{r)r^'^ = 0. (3.53) 

1 — >oo 

It follows from (j3.50p . (j3.52p . (j3.53p and integration by parts that 

1 f {p + aa) {x)L' {mp+aa (r (x ) ) )m^ (r (x) ) 



2 J r'^ix) 



1 
2 

+ 00 



+00 

a{r)L{mp+aa{r))m^{r)r~^dr 







+ / L{mpj^acr{r))m„{r)r ^dr. (3.54) 
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Since 
and 

fOO /•+00 fOO 

s^^" L{mp{s))dsa{x)dx = I a{r) i s^"^ L{mp{s))dsdr, (3.56) 

I r{x) Jo Jr 

([M]) and (f33il) imply 

"+CXD 



Q'{a) = I L{mp+aa{r))m„{r)r '^dr 

a{r) / s~^L{mp{s))dsdr. (3.57) 



Using (|3.4ip . we have m(j{r) = J^a{s)ds, so substituting this into the first term in (|3.57p 
and interchanging the order of integration gives 







+00 

L{mpj^aa{r))m^{r)r^^dr 



+00 j-T 

/ r^^L{mp+aa(.r))a{s)dsdr 
Jo 

+00 r+oo 

a{s) / L{mp+aa{r))drds 







+00 r+oo 

a{r) / s^^L{mpj^aa{s))dsdr. (3.58) 







Hence (|3371l and (|338D yield 



r-\-oo poo 

Q'{a) = / a{r) / s-3(L(mp+„<,(s)) - L{mp{s)))dsdr, (3.59) 

and therefore 

Q(0) = Q'(0) = 0. (3.60) 
Differentiating (j3.59p again, we obtain 

^ " y ^^'^^ J s"^ L' {mp+aa{s))m^{s)dsdr, (3.61) 
and interchanging the order of integration gives 

— ^-"2— = a y '^^''')'^'''^'^''^p+a(T{s))ma{s)ds. (3.62) 

Noting that a{r)dr = mo-(s), we obtain 

^^p = ay^ s'''L'{mp+Us)){m^{s)fds. (3.63) 
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Therefore, if L'{m) > for < m < M, then 

(PQ(a) 

, ; ' > 0, /or < a < 1. (3.64) 
This, together with (I5:un|) and H^Mh . yields di{p,p) = Q{1) > 0. □ 

Lemma 3.2. Let (p, v) be a solution of the Cauchy problem lll.l]) . icl. 1\) as stated in Theorem 
3.1, then 



d{p{t),p) + d^{p{t),p) - —\\V{Bp{.,t) - Bp)\\l 



E{p,^){t) - F{p) 

1 

8^ 

+ \j p{W\^ + \v'^\^){x,t)dx. (3.65) 
Proof. From Al)-A3), for any x £ Gt we have 

f{x,t)=ji{C-t{x)), (3.66) 
(see (3.26)). In view of (3.22) and (3.27), 

j\x,t)=jii^.t{x)) = L{mpMi-t{x))), (3.67) 
for x E Gt. This, together with (3.21), yields 

j2(x,t) =L(mp(t)(r(x))), x £ Gf (3.68) 
Therefore, by ()3.17p . we have 

r 1 r p(x,t)L(m„(f)(r(x)) 

E{p{t)Mt)) = J A{p){x,t)dx + -J dx 



Svr 

Here we have used the fact that 



^ j \VBp\\x,t)dx + ^ J p{\v'^\^ + \v''\^){x,t)dx. (3.69) 



/9(x,t)L(mp(i)(r(x))^_^ _ r p{x,t)L{mp^t){r{x)) 



r^{x) Jq^ r^{x) 

which holds because Dt = {x £ : p{x, t) > 0}, Gt C Dt and meas{Dt — Gt) = 0. It follows 
from (2.5) and ^^W^ that 

E{p,^){t) - F{p) 

{A{p){x,t) - A{p){x))dx 



If p{x,t)L{mp(^t){r{x)) - p{x)L{mp{r{x)) 
2 J r2(x) ^ 

-l-{\\VBp{x,t)f-\WBp\\l) 

+ \ I p{W? + \v'^?){x,t)dx. (3.70) 
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/ 



On the other hand, 

-l-{\\VBp{.Ml-\WB~p\\l) 

= -^\\V{Bp{-,t) - VB~p)\\l - ^ [ VBpix)-iVBpix,t)-VBpix))dx. (3.71) 
OTT Air J 

Noting that A{Bp — Bp) = — 47r(/> — p), and integrating by parts (this is legitimate, cf. |29j ) 
gives, 

~ i / ■ i^Bp{x,t) - VBp{x))dx 
J Bp{x){ABp{x,t) - ABp{x))dx 

Bp{x){p{x, t) - p{x))dx. (3.72) 

By (I32Q1)-(I322]), and noting (lOHl) . we have 
E{p,^)(t) - F{p) 

= I (^Aip) - Aip) + {p- p){Q ^^^^^ds - Bp}^ dx 
+ d,ip{t),p)-^m{Bpix,t)-Bp)\\l 

+ ^Jpi\v'\^ + \v^\^)ix,t)dx. (3.73) 

Since p{-,t) G Wm, f p{x,t)dx = f p{x)dx = M. Thus / X{p{x,t) — p{x))dx = 0. Therefore, 
the first term in (j3.73p is the same as d{p{t),p) defined by (3.29). This completes the proof 
of the lemma. □ 

We are now in a position to prove Theorem 3.1. 

Proof of Theorem 3.1. Assume the theorem is false. Then there exist eo > 0, t„ > and 
initial data /?„(x,0) G Wm and v„(x,0) such that for all n e N, 

d{pn{0),p) + di(po,p) + ^||Vi?p„(0) - VBpWl 

+ 11 pnixMKixM' + \vU^M'Kx)dx < (3.74) 
2 ./ n 



but for any a £ 



d{pn{tn),T''~p) + di{p{t),T''p) + ^\\VBpn{tn) " VSTVIll 

ovr 

+ \ Pn{x,tn){\vl{x,tn)\^ + \vl{x ,tn)\^){x)dx > eg. (3.75) 
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By ([3:65]) and (ITOD . we have 

lim ^(p„(0),v„(0)) =F(p). (3.76) 
Since E{pn{t),Vn{t)) is non-increasing in time, 

lim supF(p„(t„)) < lim F(p„(t„), v„(t„)) < lim F(p„(0), v„(0)) = (3.77) 

n— >oo n— >oo n— +00 

(The first inequahty holds because we have, similar to (3.71), 

Eip,v){t) - F{p{t)) = I f p{W\^ + |f3p)(x,t)dx > 0, t > 0.) 



2 „ 

Therefore {/)„(•, t„)} C Wm is a minimizing sequence for the functional F. We apply Theorem 
12.31 to conclude that there exists a sequence {a„} C M such that up to a subsequence, 

||V(5p„(t„) - Sr'^"p)||2 ^ 0, (3.78) 

as n — > oo; this is where we use the assumption that the minimizer is unique up to a vertical 
shift. Note also that for any p S Wm and a E M, 

\\VB{T-p) - VBp\\2 = \\VB{p) - VBT--p\\2, 
d{T''p,p) = d{p,T-^p), and di(Ty,p) = di{p,T-^p). 

Thus, by (j3.65p . the fact that the energy is non-increasing in time, and F{T"'p) = F{p), we 
have for any p £ Wm and a G M, 

E{pn{tn),yrn{tn))-F{T^-p) 
= d{pn {tn ) , T'^" p)+di{p{tn), T"" p) 

-l.\\V(Bpn{t^)-BT'''^p)\\l 

+ \j Pn{K\^ + \vl\'^){x,tn)dx 

<i^(/9„(o),v„(o))-F(r'^V) 

= F(/9„(0), v„(0)) - F{p) ^ 0, (3.80) 

as n ^ oo. Since 

||V5p„(t„)-V5r'^Vll2^o, 

as n ^ oo, d{pn{tn),p) > (cf. (3.31)) and di{p{tn), p) > (cf. A4) and (3.37)), we have 

d{pn{tn),T^"p) + di{p{tn),T^"p) 
+ l.\\V{Bpnitn)-T''"Bp)\\l 

+ \ I Pn{\vl,\^ + \vl\^){x,tn)dx ^ 0, (3.81) 



2 „ 

as n ^ oo. This contradicts (|3.75p . and completes the proof. 
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4 Stability of General Entropy Solutions 



In this section, we shall obtain a stability theorem for general entropy weak solutions. We 
begin with the definition of entropy weak solution. 

Definition 4.1. A weak solution (defined in Section 3) on x [0,T] is called an entropy 
weak solution of (jl.ip if it satisfies the following "entropy inequality" : 

3 3 

dtV + 5^, qj + pYj ^x,<0, (4.1) 
in the sense of distributions; i.e., 

/ / ivPt + ci-V(3-py^7]m,^^^(3] dxdt+ [ f3{x,0)ri{x,0)dx>0, (4.2) 
Jo Jr3 \ ~{ j Jr^ 

for any nonnegative Lipschitz continuous test function /3 with compact support in [0, T) x M'^. 
Here the "entropy" function t] and "entropy flux" functions qj and q, are defined by 

< q, = ^^ + m,fr-^ds=^^ + ^, (J = 1,2,3), (4.3) 

q = {qi, Q2, qs)- 

Remark 8. The inequality (j4.ip is motivated by the second law of thermodynamics ([18)). 
and plays an important role in shock wave theory ([32])- For smooth solutions, the inequality 
in (j4.ip can be replaced by equality. 

For a general entropy weak solution, our stability result is given by the following theorem: 

Theorem 4.1. Suppose 1 < 7 < 2. Let {p,m,^){x,t) (t G [0, T], x e M?) with (p, m) G 
L°°(M^ X [0, T]), be a weak solution of satisfying the entropy condition ^.1^ and let 

(p, m,l>)(x,t), t G [0,r], X eM? he any solution of satisfying (p, m) G Wl^'^i^^ x 

[0,T]). Assume 

Z{T) =: sup (||p(-,t)||oo(||p(-,t)||oo + m;t)\U^-^VolS{t)f/^ + 1 1 V,v(., t)| U) < +00, 
0<t<T 

(4.4) 

and 

— ,— G L°°(M3 X [0,T]). (4.5) 
P P 

where S{t) = Supp\p — p\{-,t). Then there is a constant C{T) depending on T and Z[T) such 
that 

Y{t) < C{T)Y{0), 0<t<T, (4.6) 
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where 



Y{t) = D{p,p){t) + ||Vp(V$ - Vl.)||i(t) + I p(x,t)|v - v|2(x,t)dx, 

$ = -Bp,^ = -Bp, 



and 



D{p, -p) 



P{p) - pip) -P'{p){p- P) 



dx. 



7-1 

Remark 9. The function (p, m, ^) in the theorem could be, but is not necessarily, a rotating 
star solution. 

Remark 10. For 1 < 7 < 2, it is easy to see 

Dip, p)>C\\p-p\\l 



for some constant C > if p G L°°(M'^ x [0, T]). 
Proof of Theorem \4-i\ 

Letting U = (p, m)"'" with m = (mi, 771,2,^3) = P^ and U = {p,m)^, we can write system 
^ as 



A<1> = Airp. 
Here the flux functions Fj[U) are given by 



(4.7) 



F2(C/) = (m2, ^,p(p) + ^, , 
Fm = {mz,^.^,p{p) + 4y . 
The entropy and entropy fluxes rj and q are as in (4.3) and satisfy 

Vqj{U) = Vr]{U)VFj{U), j = 1, 2,3, 

as is easily verifiable. Since U is an entropy weak solution 

3 3 
dtr){U) + J2 dx^qjiU) + p^, Vm, (.U)<^., < 0, 
j=i j=i 

in the sense of distributions. Because U G W^^J^ is a weak solution of (jl.ip . we have 

3 3 

j=i j=i 



(4. 



(4.9) 



(4.10) 



(4.11) 
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We define the relative entropy-entropy flux pairs by 

r]*{U,U) =r]{U) -r]{U) -Vr]{U){U -U), 

q*{U,U) = qj{U) - qj{U) -Vr]{U){Fj{U) - F,{U)) (j = 1,2,3). 
Using (|4.10p and (|4.1ip gives 

3 



(4.12) 



3 3 
= idtviU)+^d,^qj{U)) - {dt^{U) +Y,d.,q,m 

3 

3 

- V7?(i7){(C/ - U)t + Y,d.,{F,iU) - F,{U))} 

< (Vr/(C/) - Vr]{U))R - V^r]{U){R, U - U) 

3 

- V%(C7) (a,,C7, F,{U) - Fj{U) - F;{U){U - U)) , (4.13) 

in the sense of distributions, where 

R = (0, -pV^f, and R = (0, -pV^f. (4.14) 

It is easy to check that 

(V7?([/) - Vt]{U))R - V^7]{U){R, U-U) 

= -p(v-v)-(V$- V$), (4.15) 

so that 

3 

dtV*+Y.d,,qj 

< -p(v - v) • (V$ - V^) 

3 

- V^viU) ^ (a,^.C7, F,{U) - F,{U) - F;{U){U - U)) , (4.16) 
in the sense of distributions. That is, for any nonnegative, Lipschitz continuous test function 



32 



■0 on X [0, T), with compact support, we have 



> 



V'p(v - v) • (V^> - V<^)dxdt 

3 



/ / ijV^v{U)f2{dx,U, F,{U)-Fj{U)-F'{U){U -U))dxdt. 

Jo -^j^ 



A calculation gives 



and also 



So, for 1< 7 < 2, 



, p"ip) 

^ 7-1 


mi 


m2 


p^ 


mi 


1 
P 








_ m2 





1 
p 













J ) 




-p) 


1 ,^ 





7-1 



r/* > Ci(||p(-,t)|U + ||p(-,t)||oor-'(p - P)' + ^pW - > 0, 

for some positive constant ci. 

A further calculation yields, using ()4.18p . 

3 

3 

= {p{p) - pip) - P{p)iP - p)} djVj 



1 ^ 

+ - E - Vi){vj - Vj){djVi + diVj). 



Here and in the following, we use the notation: 

d= — 

Therefore, by (fil^ and (jMH), we have 

3 



j = 1, 2, 3. 



|V2ry(C7) {d^^U, F,{U) - F,{U) - F'^{U){U - U)) 
< C\\V,y{-Moorf{x,t), 
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for x£W,t£ [0,r) and some constant C > 0. Thus, KTn - K22\\ yield 

f I I dtt/jT]* + djTpq*] dxdt+ f Tp{x,0)r]* {x,0)d:i 
Jo Jr^ \ fr[ I 



> / ippiy - v) • (V^> - V^)dxdt 



ipr]*{x, t)dxdt. 




c sup iiv,.v(-,t)iu r [ 

0<t<T Jo 



(4.23) 



Using (j4.5p . it is easy to see that there exists a positive constant A, which may depend on T, 
such that 

3 

(E \Q]\')'^'{^,t) < A7i*{x,t), {x,t) e m3 X [0,T]. (4.24) 
i=i 

For fixed L > 0, t G (0, T) and small e > 0, we consider the test function ip{x, r) = (^{x, t)i9{t) 
defined by 

1, < r < t 

^i^) = { -^{t-T) + l, t<T<t + e (4.25) 
0, t + e<T <T, 



1, (x,T)eRi 
?(x,r) = <! l[L + A{t - t) - \x\] + I, (x,r)Gi?2 
0, (x,r)Gi?3, 



(4.26) 



where 



= {(x, r) : < r < T, < \x\ < L + A{t - r)}, 

i?2 = {{x, t) : < T < T, L + A{t - t) < \x\ < L + A{t - t) + e}, 

i?3 = {(x,r) :0<T<r, |x| > L + A{t-T) + e}, 

and A is the constant given in (j4.24p . Substituting this in (|4.23p . a straightforward calculation 
yields, 

1 /"^^^ /■ 

r/*(x, T)dxdT 



t J\x\<L 

< / 77*(x,0)(ix 

./|x|<L+Ai 

--ff 

e JO JL+A(t-r)<|x|<L+A(t-T)+e 
ft 



p(v - v) • (V^> - V^)dxdT 

J\x\<L+K{t-T) 

+ C sup ||V^v(-,t)||oo / / r]*{x,T)dxdT + 0{e). 

0<T<T Jo J\x\<L+A{t~T) 



(4.27) 
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The second term on the right-had side of (|4.27p is negative in view of (j4.24p . together with 
Cauchy- Schwartz inequaUty. Letting e ^ 0^ in ()4.27p gives 



/ rf{x,t)dx 

J\x\<L 

< I rj*{x,0)dx 

J\x\<L+At 
ft 



- [ [ p(v - v) • (V$ - V^)dxdT 

Jo J\x\<L+A{t-T) 

+ C sup ||V^v(.,r)||oo / / r]*{x,T)dxdT. (4.28) 

0<T<T Jo J\x\<L+A(t-T) 



We now let L +00 in (|4.28|) to get 

r]*{x, t)dx 



< J r]*{x,0)dx 

-J y p(v - v) • (V$ - V^)dxdT 

+ C sup ||V^v(-,r)||oo / [ r]*{x,T)dxdT. (4.29) 

0<r<T Jo J 

The second term on the right hand side can be estimated as follows. By Cauchy-Schwartz 
inequality, we have 

I J p{v -v) ■ {V<^ -V^){x,T)dx 

-\ j P\^~ ^)^^ + P^^*^ ~ T)dx- (4.30) 

Applying Lemma 12.31 we obtain 
J p\V<^ -V<^\'^{x,t)dx 

<\\p{.,t)\U\V{<^-^){;t\\1 

<C\\pi;T)\\^ (1 \p-p\'/^{x,t)dx^ (1 \p-p\^/^{x,T)dxy 

= C\\p{;t)\\o.( f \p-p\''^H^,t)dx] ( [ \p-p\''/Hx,r)dx] , (4.31) 



'S(r) ) \JS{r) 

where 

S(t) = supp\p- p\{-,t). 
It follows from Holder's inequality that 

p- p\^/^{x,t)dx < I [ \p- p\'^{x,T)dx] {voIS{t))^^^, (4.32) 
Sit) \Js{t) J 
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and 

\p- p\{x,T)dx] < I [ \p- p\'^{x,T)dx] {voIS{t))^/^. (4.33) 
Sit) J \Jsit) J 

Then using (|0T]) - (03]) we obtain 

p\V^-V^\\x,T)dx<C\\p{-,T)\U\\p{-,T)\\^M\pi^ 

(4.34) 

In view of (g^Ql), (4.29), (4.30) and (4.34), we have 

J r]*{x,t)dx< J ri*{x,0)dx + CZ{T) J rf{x,T)dxdT, (4.35) 
for < t < T, where 

Z{T) =: supo<t<T{\\p{-Moo{M-Moo + I !/?(•, t)||oo)'"^(V^o/S(t))2/3 + ||V,v(-,t)||oo). 

Then (4.6) follows from Gronwall's inequality applied to (4.35) and using (4.19) and (4.20). 
This completes the proof of Theorem 4.1. 

5 Uniform A Priori Estimates 

The theorem proved in this section gives a uniform a priori estimate for the entropy weak 
solution defined in (4.2) of the Cauchy problem p.l|) and (|3.ip . As we shall see, this estimate 
justifies some assumptions made in Section 3 and should be useful for obtaining the existence 
of global weak solutions for the Cauchy problem. 

Theorem 5.1. If{p,m) G L°°([0, T]; L^(M'^)) satisfies the first equation in (1.1) in the sense 
of distributions, then 

/ p{x,t)dx= / p{x,0)dx =: M, 0<t<T. (5.1) 

Let (p, m, <5) be a weak solution defined in Definition 3.1. Suppose (/>, m, <I>) satisfies the 
entropy condition (4.2), p G L°°{[0,T]; L^{R^)) n L'^{[0,T]; U' {R^)) for some r satisfying 
r > 3/2 andr>j,m£ L°°([0, T]; L*(IR3)) (s > 3^, (r?,q) G L°°([0, T]; (M^)), where rj and 
q are given in (4.3). Moreover, we assume that (/?, m) has the following additional regularity: 

lim / / \p{x,T + h) - p{x,T)\dxdT = 0, te{0,T),a.e. (5.2) 

Then 

E{t) < E{0), 0<t<T, (5.3) 
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and if "f > |, then 

H{t) < CiH{0) + C2, 0<t<T, (5.4) 
where Ci and C2 are two positive constants only depending on 7 and AI ( cf. h5. 1\) ), where 

Jr3 7-1 Ip 87r 



Remark 11. (jS.ip and (j5.3p justify some assumptions made in Section 3 on the conservation 
of total mass and non-increase of energy. 

Remark 12. The boundedness of p'^ {x,t)dx was proved in [lOj for smooth solutions if 
7 > 4/3. Here we prove that this is stih true for general week solutions satisfying the entropy 
condition even without assuming that p G In fact, the global existence of radial L°°- 

solutions was proved in [34J for (jl.ip outside a ball. The blow-up of L°°-norm of the radial 
solutions of (jl.ip in the entire space was discussed in [27] and [11], respectively. 

Remark 13. Condition (|5.2p can be assured by the following condition 

lim sup / \p{x,t) — p{x — ey,T)dx = 0, (5-5) 

if (p,m) G L°°([0,r];L^(R3)); this is proved in the Appendix. Note that ([52]) is the 
modulus of continuity in time and ()5.5p is the modulus of continuity in space. 

In order to prove this theorem, we begin with the following lemma. 

(r > 1), then 




with 1/p = 1/r - 2/3, if r < 3/2, 
if r > 3/2; 



with 1/g = 1/r - 1/3, if r < 3, 
if r > 3. 



(5.6) 



(5.7) 



The proof of this lemma follows from the extended Young inequality (cf. [28], p. 32). 
Lemma 5.2. Suppose < pG L°°{[0,T]; L'^{M.^)) and G L°°{[0,T];L'^{M.^)), then 

m G L°°([0,r];Li(M3)). (5.8) 
Proof. Using Holder inequality, we have 

J \ui\dx = J W-^dx < ( I pdxfl\ I ^)i/2. (5.9) 
Note that ([SJD implies m G L°°([0, T]; ^^(M^)). □ 
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Remark 14. r? G L°°([0, T]; L1(M3)) implies ^ G L°°([0, T]; L^l 

Lemma 5.3. Let (p, m, <1>) 6e a weak solution defined in Definition 3.1. Suppose (p, m, <I>) 
satisfies the entropy condition (4.2), p G L'=^{[0,T]; L^{R^)) n L°^([0, T]; ^''(M^)) /or some r 
satisfying r > 3/2 and r > j, m e {[0,T]; L' (R^)) (s > 3), (r?,q) G L°°([0,r];L^ 
where rj and q are given in (4-3). Then, for any r G [0, T), we /laue 



ri{x,T)dx — / / m • V<I>dx(ii < / r]{x,0)dx, rG(0,T),a.e. 
Proof. For a fixed r G (0, T), and small positive e and ii > 0, we define 



(5.10) 



9{t) 



< t < r. 



-i(t-r) + l, r<t<T + e, 



(5.11) 



and for x G M'^, 



q(x) = < 



T + e<t<T, 



\x\ < R, 



-^{\x\-R) + l, R<\x\<R + e, 



(5.12) 



0, 



Ixl > R + e. 



Let (3{x,t) = 9{t)a{x), then /3(x,t) is Lipschitz continuous, with compact support in [0,T) x 
R?. Using (|4.2p . a calculation yields 



/■T+e 



f r]{x,t)dxdt / / r]{x,t)a{x)dxdt 

\x\<R e JR<\x\<R+e 

3 

(^gj-^)6'(t)dxdt 



+ / ri{x,0)dx + / 7y(j;, 0)a(x)(ix 

/|a;|<_R JR<\x\<R+t 



nrr+e r- 
m-V^dxdt+ / m-\/^f3{x,t)dxdt>0. (5.13) 

,-c|</? Jr J|a:|<i?,+e 



Since (r?,q) G L°°([0, T], L1(R3), we have 



lim / r]{x,t)a{x)dx = 0, a.e., tG[0, T], 

R^oo Jfi<\x\<R+e 

lim / r]{x,0)a{x)dx = 0, 

R^'^ J R<\x\<R+€ 



and 



lim 



(^gjpr)6'(t)d3; = 0, a.e., t G [0,r]. 



(5.14) 



(5.15) 



(5.16) 
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We let ^ oo in (l5l^ to get 
1 



e 



r]{x,t)dxdt + I ri{x,0)dx 



+ / / m-\/^dxdt+ / / m-V^Pix,t)dxdt>0. (5.17) 
Jo Jm.^ Jt Jr^ 



Because p G L°°{[0,T]; U (R-^)) with r > 3/2, by (5.7) we have V$ G L°°([0, T]; L''(M'^) with 
g > 3. It then fohows from (5.8), the assumption m e L°°([0, T]; L'*(M^)) with s > 3 and 
Holder's inequality that 

m- V«> G L°°{[0,T];L^{E.^)). 

This implies 

lim / / m • V$/3(x, t)(ix(it = 0. 

Letting e ^ in (ISTTj) . we obtain (l5l^ . □ 

Lemma 5.4. Lei (p, m, $) 6e an entropy weak solution defined in Section 4 satisfying the 
conditions in Lemma 5.3. Then 

dMx, t) = - [ m(y, t) ■ Vy{-^)dy. (5.18) 

Jr3 \y-x\ 

Moreover 

dt^ £ L'^{[0,T];L'^{R^)), (5.19) 

and 

G L°°([0,r] X M^). (5.20) 

Proof. The key is to prove (j5.18p . Once (|5.18p is proved, (j5.19p and (|5.20p follow from the 
fact that m G L°°{[0,T]; L'^{R^)) n L°°([0, T]; L^(M3)) and the extended Young's inequality 
(cf. [28], p. 32). In order to prove (j5.18p . we use the fact that (p, m) satisfies the first equation 
of (jl.ip in the sense of distributions. For this purpose, we choose a C°° function 6{z) {z G M^) 
with compact support in the interval [1,2] satisfying < 6{z) < 1 and 6{z)dz = 1, and 
let 

5^{z) = -<5(-), a,{z) = I 6,{s)ds, z G R^, (5.21) 

for small positive e. For y G M'^, < e < ^ and i? > 1, we set 

n I aMy\), |y| < # + e, 

[ a,{R + 2€-\y\), |y|>|+e. 

Then 

fR[y) = 0, as \y\ < e, or \y\> R + e, 

0</f(y)<l, as e < |y| <2e, or ii< |y| <i? + e, (5.23) 
/,^(y) = 1, as 2e < |y| < R. 
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d_ 
di 



For X E M^ we choose 

g^{y) = f,%-x)-^. (5.24) 

\y - x\ 

Then g^{y) G Cq"(M^) for any fixed x G M'^. Since (/f, m) satisfies the first equation of (jl.ip 
in the sense of distributions, it is easy to show (see [E] for instance), f^s p{y,t)g^{y)dy is 
differentiable in t for t G [0,T] a. e., and satisfies 

/ p{y,t)g^{y)dy= [ m{y,t) ■Vygf{y)dy, te[0,r], a. e. (5.25) 

We also let 

<7.(y) = lim 5f (y), ye Ml (526) 

it— >oo 

Then we show (|5.18p in the following steps. 

Step 1. We show that jj^g p{y,t)g^{y)dy is differentiable for t £ (0, T], a.e., and 

4/ P{y,t)ge{y)dy = I in-Vyg^{y)dy, (5.27) 

for t G (0,T],a.e. 

For this purpose, we prove that 

1 f ,(,,t + /.)-,(,,t) ^,^^^^^^ f 

'1 Jr3 ft J^3 

as /i ^ uniformly in i? for R> 1. (5.28) 

This is proved as follows. Since {p, m) satisfies the first equation of (jl.ip in the sense of 
distributions and g^{y) G C^(IR^), it is easy to verify (see |15) for instance), 

/ ipiy,t + h)- piy,t))g^{y)dy= [ [ ni{y,s)-Vyg^{y)dyds, (5.29) 
for [t,t + h] C [0,r]. Thus, 

1 f iP^yil±R^P^gn^y)dy = hm i r' [ m(y, .) . Vyg^iy)dyds. (5.30) 
/i^o JiR3 ft h^o ft 7i 7^3 

On the other hand, 



■t+h 



1 



^{y, s) ■ Vyg^{y)dyds - / ni{y, t) ■ Vygf{y)dy 
(m(?/,s) - m(y,t)) • Vygf{y)dyds 



m( 

i / (m(2/, s) - m(y, t)) • V,( "f ^ ~ f )rf;/d5 

"-A Je<\y-x\<R \y-x\ 
rt+h 



+ \ f [ (m(y, s) - m(y, t)) • Vyg^{y)dyds. (5.31) 

-J R<\y-x\<R+€ 
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The first term can be handled as follows. For h > 0, 

r-t+h 



\\l I (m(y,s)-m(y,t)).V,(^f^^)dyds| 

hJt Je<\y-x\<R \y-x\ 

h Jt Je<\v-x\<R \y-x\ \y-x\^ 



2 f 1 

< — - / / \m{y,s) -m{y,t)\- -dyds 

Jt Je<\y-x\<2e W - x\ 

1 f 1 

+ -/ / \m{y,s) -m{y,t)\- r^dyds 

h Jt J2e<\y-x\<R \y - x]'' 

2 f 

<^/ / \ui{y,s) -ui{y,t)\dyds 

^ 1^ Jt Je<\y~x\<2e 
^ i-t+h j- 

+ TU: / \m{y,s) -ui{y,t)\dyds. (5.32) 

4e n Jt J2e<\v-x\<R 



The last term in (j5.3ip can be estimated as follows. 



t+h 

(m(y,s) - m{y,t)) ■ Vyg^{y)dyds\ 

R<\y-x\<R+e 



1 11 /■*+^ f 1 

<(- + ^)t/ / \m{y,s) -m{y,t)\- -dyds 

e R h Jt jR<\y_^\<R+^ \y - x\ 

+ / \miy,s)-iniy,t)\dyds. (5.33) 

e ^ n Jt Jm.^ 

Since we choose R>1, dOT]) . ([02]) and ([03]) yield, 

/ /" m{y,s)-Vyg^{y)dyds- f m{y,t) ■Vygf{y)dy\ 

+ - + /" |m(y,s) -m(2/,t)|d2/(is. (5.34) 



Since m G L°°([0, T]; L1(M3), we have 

rt+h 

iim — I 

h 



1 /■*+^ /■ 

lim - / / \m{y,s) -m{y,t)\dyds = 0, t e [0,T], a.e. (5.35) 
.^0+ Jt Jks 



Therefore ^ //^^ /ig3(m(?/, s) — m(y, t)) • V ygf{y)dyds converges to zero as h ^ 0+ uniformly 
in it! for i? > 1. By a similar approach, we can show that ^ fll_^ f^s{m{y, s) — m(y,t)) • 
\7yg^{y)dyds converges to zero as /i — > 0— uniformly in R for > 1. This verifies (j5.28p . 
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(j5.27p follows by the following argument, using (5.22) and (5.24). 

4 / P{v,i)9e{v)dy 

= lim lim / ^^y^^±4^^g^iy)dy 

= lim lim / ^^^^^^±4^^g^{y)dy 
d C 

= i™ 17 / piy^'t^)9^iy)dy 

it— >oo at _/jg3 
= lim / m • Vygf{y)dy 

™.-Vyge{y)dy. (5.36) 

Step 2. In this step, we show that 



p{y,t)ge{y)dy ^ [ f^^^dy as e ^ 

Jr3 \y - x\ 

uniformly in t for t G (0, T), (5.37) 

and 

linr / m{y,t) ■Vyge{y)dy = m{y,t)-Vy{- r)dy, (5.38) 

for t G (0,r). 

We prove (lOTD as follows. Since p G L°°([0, T]; L'^(M3)) with r > 3/2 and r > 7, we have, 
by using Holder inequality, 

'^'^'\dy- [ piy,t)gMdy\< f ^dy 



< 



y-X\ J]R3 J e<\y-x\<2e \y - X\ 

p^{y,t)dyfl'{f V^M^' 

<^<\y-x\<2t Jt<\y~x\<2t\y ~ X\ 



.2e l/« 



< IIOI|ioo([o,T];L' 



)( / 47rs2~'ds) , (5.39) 
where I = Since r > 3/2, Z < 3, (|5.37p follows. Next (|5.38p can be shown as follows. 
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Since m E L°°{[0,T];L'{R^)) for s > 3, we have 

1/ m{y,t)-\7yg^{y)dy- I m{y,t)-\7y- — - — rdy\ 

jR3 \y-x\ 

= \ • Vy(5(,(y) - I 

Je<\y-x\<2e lU — ^\ 

= 1/ in{y,t)-Vy{-^ia,{\y-x\)-l))dy\ 

Je<\v-x\<2e W ~ ^1 



< / \ni{y,t)\i-^— + -^5,i\y-x\)))dy 

Je<\y~x\<2e lU ~ ^\ W ~ ^\ 

2 f 1 

< - / I rdy 

£ Je<\y-x\<2e W ~ ^1 



2,, 

< — m 

e 



Je<|j/-x|<2e |y ~ ^1 

< -||m||Loo([o,r];L'>(R3)) ( / 47rT^"^'dT , (5.40) 



where q = Since s > 3, then q < 3/2. Therefore, (|5.38p is proved. 

By ()5.37p and (5.24), we have that J^^ j^z^dy is differentiable with respect to t for {t,x) £ 
(0,r) X a. e. Moreover, by (5.24), (|07D and dOH]) . we obtain, 

^/ r^dy = ^(lim [ p{y,t)gMdy) 

= f P(y ^ 9e{y)dy = lim [ m{y,t) ■ Vyge{y)dy 

= [ in{y,t)-Vy{-^)dy. (5.41) 

Jrs \y-x\ 

This proves (5.18). (5.19) and (5.20) then follows as we showed at the beginning of the proof 
of this Lemma. □ 

Proof of Theorem 5. 1 
We prove Theorem 5.1 in the following steps. 

Step 1 In this step, we prove (5.1). This can be proved by using (5.25) in which gf'{y) is 
replaced by (y), i.e., 

^/ p{y,t)f^{y)dy= f in{y,t) ■ Vyf,''{y)dy,t e [0,T],a.e., (5.42) 
where f^ is defined in (5.22). We integrate (|5.42p to get 

/ p{y,t)f^{y)dy- f p{y,0)f^{y)dy = f f ui{y, s) ■ Vyf^{y)dy, (5.43) 
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By using a same argument as in the proof of Lemma 5.4, we can prove 

lini lim / p{y,t)f^{y)dy = p{y,t)dy, 



lini lim / piy,0)f^{y)dy = j p{y,0)dy, 



and 



lini lim / / m{y, s) ■ Vyf^{y)dy = 0. 
(5.1) follows from (j5.43p by letting R ^ oo and e ^ 0. 
Step 2 In this step, we show that 



/ / p{x,s)ds^{x,s)dxds = ]-( ! {p<^){x,t)dx- [ {p^){x,0)dx\ t e [0,T). (5.44) 
Jo Jr^ ^ \Jr3 Jr3 J 



This is can be proved as follows. 




p{x^ s)ds^{x, s)dxds 




JR3 

lim — 

h^o h Jq 

lim — 

h^O h \J;^ 

lim — 

/i-+o h Ji^ 



p{x,s) 



t+h 




p{y,s + h) - piy,s) 
\x - y\ 
p{x,s- h)p{y,s) 



dydxds 




dydxds 

x-y\ 

{p{x,s- h) - p{x,s))p{y,s) 




p(x.)p(,,.) ^ \ 
\x-y\ ) 



lim — 

/i^o h 

■ lim — 

h^O h 



t+h 



\x - y\ 
p{x,s- h)p{y,s) 



dydxds 



1 




\x - y\ 

{p{x,s- h) - p{x,s))p{y,s) 



dydxds — lim , , 

/l-»0 h Jq 



p{x,s - h)p{y,s) 
\x - y\ 



\x - y\ 

+ / {p^){x,t)dx- I {p<^){x,0)dx. 



dydxds 



dydxds 



(5.45) 



On the other hand. 



1 



h Jh 7k3 
_ 1 
~ h 
1 



{p{x,s- h) - p{x,s))p{y,s) 
\x - y\ 

(p(x,r) - p{x,T + h))p{y,T) 



dydxds 
dydxdr 



+ 



10 



t 



\x - y\ 

{p{x, r) - p{x, T + h)){p{y, T + h)- p{y, r)) 
\x - y\ 

{p{x, t) - p{x, T + h))p{y, T + h) 



dydxdT 



Since 



lim — 

h^O h 



\x -y\ 

1 f p{y,T + h) - p{y,T) 



\x - y\ 
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dydxdT. 



dy = —dr^{x, t) 



(5.46) 



1 1 r P(y,r^k) p(y,r) ^^^ ^ ^ ^ ^^^^^^ 

IT- Jr3 \x - y\ 

for small \h\. Therefore, 

,1 /■* /■ f {p{x,t)- p{x,T + h)){p{y,T + h)- p{y,T)) 



h Jo 7]R3 7iR3 \x - y\ 

< (Pt^llL-([o,T]xR3 + 1) / / \{p{y,T + h) - p{x,T)\dydT. 

Jo Jr3 



dydxdT\ 



Then ([O]), (5.20) and (|5:i8]) imply 



Similarly, we have, for small 

1 /■* f f {p{x,t) - p{x,T + h))p{y,T + h) 



t-h Jr3 f ~ y| 



dydxdrl 



(5.48) 



'^--o /i Jo Jr^ Jr3 f - y| 



< (Pt^llL-([o,T]xR3 + 1) / [ p{y,T + h)dydT. (5.50) 
Since p G L°°([0, T]; L^M^))^ ([530D implies, 

h~*o h Jt^hJR^ Jr3 \x-y\ 
Hence, (fSliGl) . (I5:i9l) and ([53T|) yield 

limi/"7 / ^^^^^^^4^#^^^^d2/dxcir 

= hm i /7 / ^"^"'"^""^"'"|^^^"^"'"^ ci,dxdr. (5.52) 
/^-^o Ai Jq J^3 J^s \x - y\ 

This, together with (|5.45p . imphes (|5.44p . 

Step 3 In this step, we prove (5.3). 

Since p £ L°°{[0,T; L^{R^)) f] L°°{[0,T]; L"" (R^)), where r > 3/2 and r > 7, we have, in view 
of (5.7) that 

V$ e L°^{[0,T];L^/^{R^))nL'=°{[0,T];L^{W^)), (5.53) 
if r < 3, where j = ^ ~ We also know that A > 3 if r > 3/2. Similarly, by (5.7), we have 

V$ G L°°{[0,T];L^/^{R^))nL'^{[0,T] x R^), (5.54) 

if r > 3. Furthermore, because (yO, m) satisfies the first equation of (1.1) in the sense of 
distributions, then by a density argument as in [15], in view of (5.19), (5.20), (j5.53l and 
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(|5.54p . we have, 

= / p{x, s)ds^{x, s)dxds + / / m{x, s) ■ ^/^{x, s)dxds, (5.55) 
Jo Jo Jr^ 

for t £ [0, T). This, together with (5.10) and (5.44), imphes (5.3), due to the fact 

E{t) = I r]{x,t)dx-- I {p<!>){x,t)dx, (5.56) 



2 

for t £ [0,T). 

Step 4 In this step, we proof (5.4). 

First, since p G L°°{[0,T]; L^{R^) D L°°{[0,T]; L'' (M.^) with r > 3/2, it follows from 
and [30] that 

If ip^){x,t)dx = -^ [ \V^\x,t)dx, t£[0,T]. 
Using (2.19), we have, for 7 > 4/3 

^\V<l>\'^}dx = j ^pBpdx <C j p'^'^dxi j pdxf'^ = M^/^ j p'^'^dx, (5.57) 

-7-1 

where A{p) is given by (2.3). Taking p = 1, q = 4/3, r = 7, and a = hi Young's 

inequality (2.17), we obtain, 

IHl4/3<IHI?llp||^-'^ = Ml|p||^-^ (5.58) 

This is 

J p^/^dx < Mh{ J p^dx)\ (5.59) 
where h = 3(:^~ij- Since 7 > 4/3, we have < 6 < 1. Therefore, (15.570 and ()5.59p imply 

j ^pBpdx < C(7 - lfMh+l{ J A{p)dx)K (5.60) 
Using the inequality (cf.[15] p. 145) 

aP < ea' + e"*/'/?*, (5.61) 

if + = 1 (s, f > 1) and e > 0, since 6 < 1, we can bound C{j - l)^Mt"+i (/ A{p)dxf 
by ^ / A{p)dx-'rC2, where C2 is a constant depending only on M and 7 (we can take e = 1/2 
and s = 1/6 and t = {I — s^^)^^ in (I2.26P since s > 1 due to < 6 < 1). Therefore, 

h [ {p^){x,t)dx\ = ^\\VH;t)\\l<l [ ^^^dx + C, (5.62) 
2 Jr3 87r 2 J]g3 7 - 1 

for t G [0,T), where C is a constant only depending on M = J^^ j p{x,t)dx = J^^ p{x,0)dx 
(cf. (5.1)) and 7. This, together with (5.3), implies (5.4). 
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6 Appendix 

In this appendix, we prove the following theorem which is Remark 1131 in Section 5. 



Theorem If (/O, m) G L°°([0, T]; L^(M'^)) satisfies the first equation of (jl.ip in the sense of 
distributions, then 



\\y\ 

implies 



lim sup / \p{x,t) — p{x — ey,t)\dx = ^, t € (0, T), o.e., 
^^Oo<t<T,|i;|<i Jr3 



(6.1) 



lim / / \p{x,t + h) - p{x,t)\dx = {), te{0,T),a.e. (6.2) 
Proof. For any fixed t G (0, T) and small h, we let 

w{x) = p{x, t + h) — p{x, t). 
First, we note that if ip{x) £ C^(M^)) with ip and Vip being bounded in M'^, then 

rt+h r- 

w{x)'il>{x)dx = I I m{x, s) ■ Vip{x)dxds. (6.3) 



This is because {p, m) G L°°([0, T]; L^(M^)) satisfies the first equation of (jl.ip in the sense of 
distributions. The justification of (j6.3p is standard, for instance, see pl5J. In view of (|6.3p . 
we have 

w{x)'il^{x)dx\ < h sup |VV'(x)|||m||i.o([o,r];Li(R3))- (6.4) 



We choose as 



ipix) = / sgn{x - ey)5{y)dy, 



where sgn is the sign function, 5 G C^(M^) is a smooth function satisfying < 5{y) < 
1,J^3 6{y)dy = 1 and supp S C {y G : \y\ < 1}. Then < ^ for some constant C. 

Moreover, 

I / w{x)-il;{x)dx — / |u)(x)|(ix| 
Jr« Jrs 

{w{x) — w{x — ey))sgn{x — ey)5{y)dydx\ 



< sup / \w{x) — w{x — ey)\dy 

\y\<l Jr3 

< sup / \p{x,t) — p{x — ey,t)\dx + sup / \p{x,t + h) — p{x — ey,t + h)\dx. (6.5) 

\y\<lJR3 \y\<lJR3 
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Therefore, 

\'w{x)\dx 



< sup / \p{x,t) — p{x — ey,t)\dx + sup / \p{x,t + h) — p{x — ey,t + h)\dx 
|j/|<iJr3 \y\<iJK^ 

Ch 

+ — ll"i||L°°([0,T];Li(M3))- (O-O) 

We let /i ^ first in (j6.6p . (j6.2p follows from (j6.1|) because e is arbitrary. □ 
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